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We generalize the typical medium dynamical cluster approximation to multiband disordered sys¬ 
tems. Using our extended formalism, we perform a systematic study of the non-local correlation 
effects induced by disorder on the density of states and the mobility edge of the three-dimensional 
two-band Anderson model. We include inter-band and intra-band hopping and an intra-band dis¬ 
order potential. Our results are consistent with the ones obtained by the transfer matrix and the 
kernel polynomial methods. We apply the method to Ka;Fe 2 -;/Se 2 with Fe vacancies. Despite the 


strong vacancy disorder and anisotropy, we find 
results demonstrate the application of the typical 
to study Anderson localization in real materials. 

PACS numbers: 71.23.An,72.80.Ng,71.10.Fd,74.70.-b 

I. INTRODUCTION 

The role of disorder (randomness) in materials has 
been at the forefront of current researcbi'^— due to the 
new and improved functionalities that can be achieved 
in materials by carefully controlling the concentration of 
impurities in the host. At half-filling and in the absence 
of any spontaneous symmetry breaking field, disorder 
can induce a transition in a non-degenerate electronic 
three-dimensional system from a metal to an insulator 
(MIT)ii^. This phenomenon, which occurs due to the 
multiple scattering of charge carriers off random impuri¬ 
ties, is known as Anderson localization^. 

The most commonly used mean-field theory to study 
disordered systems is the coherent potential approxima¬ 
tion (CPA)^“— , which maps the original disordered lattice 
to an impurity embedded in an effective medium. The 
CPA successfully describes some one-particle properties, 
such as the average density of states (ADOS) in substi¬ 
tutional disordered alloys^^— . However, being a single¬ 
site approximation, the CPA by construction neglects all 
disorder-induced nonlocal correlations involving multiple 
scattering processes. To remedy this, cluster extensions 
of the CPA such as the dynamical cluster approximation 
(DCA)S“— and the molecular CPAji^ have been devel¬ 
oped, where nonlocal effects are incorporated. Unfortu¬ 
nately, all of these methods fail to capture the Anderson 
localization transition since the ADOS utilized in these 
approaches is neither critical at the transition or distin¬ 
guish the extended and the localized states. 

In order to describe the Anderson transition in such 
effective medium theories, a proper order parameter has 
to be used. As noted by Anderson, the probability dis- 
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tribution of the local density of states (LDOS) must be 
considered, and the most probable or typical value would 
characterize it^d^. It was found that the geometric mean 
of the LDOS is a good approximation of its typical value 
(TDOS) and it is critical at the transitioni^^— , which 
makes it an appropriate order parameter to describe An¬ 
derson localization. Based on this idea, Dobrosavljevic 
et ali^ formulated a single-site typical medium theory 
(TMT) for Anderson localization which gives a qualita¬ 
tive description of the transition in three dimensions. In 
contrast to the CPA, the TMT uses the geometrical av¬ 
eraging over the disorder configuration in the self consis¬ 
tency loop. And thus, the typical not the average DOS is 
used as the order parameter. However, due to the single¬ 
site nature of the TMT it neglects nonlocal correlations 
such as the effect of coherent back scattering. Thus, the 
TMT underestimates the critical disorder strength of the 
Anderson localization transition and fails to capture the 
reentrant behavior of the mobility edge (which separates 
the extended and localized states) for uniform box disor¬ 
der. 

Recently, a cluster extension of TMT was developed, 
named the typical medium dynamical cluster approxima¬ 
tion (TMDCA)i^, which predicts accurate critical disor¬ 
der strengths and captures the reentrant behavior of the 
mobility edge. The TMDCA was also extended to include 
off-diagonal in addition to diagonal disorderi^S. How¬ 
ever, like the TMT, the previous TMDCA implementa¬ 
tions have only been developed for single-band systems, 
and in real materials, there are usually more than one 
band close to the Fermi level. Sen performed CPA calcu¬ 
lation on two-band semiconducting binary alloys^S, and 
the electronic structure of disordered systems with multi- 
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pie bands has also been studied numerically in finite sys- 
tems2ii2^. But a good effective medium theory to study 
Anderson localization transition in multiband systems is 
still needed to understand the localization phenomenon 
in real systems such as diluted doped semi-conductors, 
disordered systems with strong spin-orbital coupling, etc. 

In this paper, we extend the TMDCA to multiple 
band disordered systems with both intra-band and inter¬ 
band hopping, and study the effect of intra-band disor¬ 
der potential on electron localization. We perform cal¬ 
culations for both single-site and finite size clusters, and 
compare the results with those from numerically exact 
methods, including transfer matrix method (TMM) and 
kernel polynomial method (KPM). We show that finite 
sized clusters are necessary to include the nonlocal ef¬ 
fects and produce more accurate results. Since these re¬ 
sults show that the method is accurate and systematic, 
we then apply it to study the iron selenide superconduc¬ 
tor Ka;Fe 2 -ySe 2 with Fe vacancies, as an example to show 
that this method can be used to study localization effects 
in real materials. In addition, as an effective medium the¬ 
ory, our method is also able to treat interactions^^, unlike 
the TMM and KPM. 

The paper is organized as follows. We present the 
model and describe the details of the formalism in Sec. m 
In Sec. IIII A[ we present our results of the ADOS and 
TDOS for a two-band disordered system with various 
parameters, and use the vanishing of the TDOS to: de¬ 
termine the critical disorder strength, extract the mo¬ 
bility edge and construct a complete phase diagram in 
the disorder-energy parameter space for different inter¬ 
band hopping. In Sec. ME we discuss simulations of 
K 2 ;Fe 2 -ySe 2 with Fe vacancies. We summarize and dis¬ 
cuss future directions in Sec. CYl In Appendix El we 
provide justification for the use of our order parameter 
ansatz. 


II. FORMALISM 

A. Dynamical cluster approximation for multiband 
disordered systems 

We consider the multiband Anderson model of non¬ 
interacting electrons with nearest neighbor hopping and 
random on-site potentials. The Hamiltonian is given by 

<ij> a,P=l 

+E E 

i=l a,/3—1 

The first term provides a realistic multiband description 
of the host valence bands. The labels i, j are site indices 
and a, P are band indices. The operators cl^icia) create 
(annihilate) a quasiparticle on site i and band a. The 


second part denotes the disorder, which is modeled by a 
local potential that is randomly distributed accord¬ 
ing to some specified probability distribution 
where fi is the chemical potential, and 

are the hopping matrix elements. Here we consider bi¬ 
nary disorder, where the random on-site potentials 
obey independent binary probability distribution func¬ 
tions with the form 

p{vr^)=xsivr^ - vf )+(1 - x)5ivr^ - vf). ( 2 ) 


In our model, there are lb band indices so that both 
the hopping and disorder potential are 4 x lb matrices. 
The random potential is 


/ Faa . . . y o^P \ 


Yi = 

\vp 

while the hopping matrix is 

/ /Q'Q: 

tij = 


\*^J 


vfV 




fia ,PP 




( 3 ) 


( 4 ) 


where underbar denotes It x lb matrix, and are 
intra-band hoppings, while and are inter-band 
hoppings. Similar definitions apply to the disorder po¬ 
tentials. If we restrict the matrix elements to be real, 
Hermiticity requires both matrices to be symmetric, i.e., 
r/3 = and 

To solve the Hamiltonian of Eq. (TJ we first generalize 
the standard DCA to a multiband system. Within DCA 
the original lattice model is mapped onto a cluster of size 
Nc = with periodic boundary condition embedded in 
an effective medium. The first Brillouin zone is divided 
in Nc coarse grained cellsi^, whose center is labeled by 
K, surrounded by points labeled by k within the cell. 
Therefore, all the A:-points are expressed as k = K + k. 
The effective medium is characterized by the hybridiza¬ 
tion function A{K, w). The generalization of the DCA to 
a multiband system entails representing all the quantities 
in momentum space as 4 x lb matrices. 

The DCA self-consistency loop starts with an initial 
guess for the hybridization matrix A(A', w), which is 
given by 


A(/C,w) = 


/A““(A:,u;) ••• A°P{K,oj)\ 


\A'5“(K,a;) ••• APP{K,uj) ) 


( 5 ) 


For the disordered system, we must solve the cluster 
problem in real space. In that regard, for each disorder 
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configuration described by the disorder potential V we 
calculate the corresponding cluster Green function which 
is now an IbJVc x Ib-^c matrix 

G^(V) = (wl - . (6) 

Here, I is identity matrix and A^- is the Fourier transform 
(FT) of the hybridization, i.e., 

A°‘^(K)exp[iK ■ (r, - r^)]. (7) 

K 


Green function. 

The DGA self-consistency condition requires the dis¬ 
order averaged cluster Green function equal the coarse 
grained lattice Green function 

Gc{K,u;) = GiK,u;) . (12) 

Then, we close our self-consistency loop by updating 
the hybridization function matrix using linear mixing 

A„(jF,c^) = A,{K,uj)+^[G:\K,uj)-G-\K,w)], (13) 


We then stochastically sample random configurations 
of the disorder potential V and average over disorder 
((• • •)) to get the IbNc x IbNc disorder averaged cluster 
Green function in real space 


/(Gr(w,^))y ••• {Gf{u:,V))^^\ 




V(Gf“(^,H)),^. ••• {GP/{u^,V))^J 

( 8 ) 

We then Fourier transform to K space and also impose 
translational symmetry to construct the 7C-dependent 
disorder averaged cluster Green function Gc{K, w), which 
is a 4 X lb matrix for each K component 


Gc{K,lo) = 


/GT{K,uj) ••• Gf[K,u:)\ 


VGf“(iF,w) ••• GP/{K,u)) 


(9) 


After the cluster problem is solved, we can calculate the 
coarse grained lattice Green function matrix 

/g““(a:,w) ••• g“^(a:,w)\ 


where the subscript “o” and “n” denote old and new re¬ 
spectively, and is a linear mixing factor 0 < .^ < 1. 
The procedure above is repeated until the hybridiza¬ 
tion function matrix converges to the desirable accuracy 
Ar,{K,uj) = A,{K,lo) . 

We can see that when the inter-band hopping, , and 
disorder potential, 17“^, vanish all the h x 4 matrices 
become diagonal, and the formalism reduces to single 
band DGA for lb independent bands. 


B. Typical medium theory for multiband 
disordered systems 

To study localization in multiband systems, we gener¬ 
alize the recently developed TMDCAi^ where the TDOS 
is used as the order parameter of the Anderson localiza¬ 
tion transition, so the electron localization is captured 
by the vanishing of the TDOS. We will use this TMDCA 
formalism to address the question of localization and mo¬ 
bility edge evolution in the multiband model. 

Unlike the standard DGA, where the Green function is 
averaged over disorder algebraically, the TMDGA calcu¬ 
lates the typical (geometrically) averaged cluster density 
of states in the self-consistency loop as 


G{K,u:) = 


(10) . V (14) 


\G^"iK,uj) ••• G^^{K,w)J 
= ^ ( GcjK, + A{K, uj) -Sk + e(K)) ^ 

k 

where the overbar denotes cluster coarse-graining, and 
e{K) is the cluster coarse-graining Fourier transform of 
the kinetic energy 

+ ( 11 ) 

k 

where Eq^ is a local energy, which is used to shift the 
bands. The diagonal components of Eq. [TU] have the 
same normalization than a conventional, i.e., scalar. 


which is constructed as a product of the geometric av¬ 
erage of the local density of states, pu = —^ImGu^uj), 
and the linear average of the normalized momentum re¬ 
solved density of states p{K,uj) = —^ImGc{K,uj). The 
cluster-averaged typical Green function is constructed via 
the Hilbert transformation 


GlyP{K,uj)= j 


doj 


OJ — oj' 


(15) 


Generalization of the TMDGA to the multiband case 
is not straightforward since the off-diagonal LDOS 

>^01 positive definite. We con¬ 
struct the 4 X lb matrix for the typical density of states 
as 
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Ptyp 


eTv^ 


{u))) 


T^E.prM 




t^E Jpr(<^ 


eT^E.(Hp“/Ml) / p°'^(^.^) \ \ 




(16) 


The diagonal part takes the same form as the single-band 
TMDCA ansatz, and the off-diagonal part takes a similar 
form but involves the absolute value of the off-diagonal 
‘local’ density of states. 

We construct the typical cluster Green function 
through a Hilbert transformation 




GlJK,uj) = 




Jduj 


/ Ptl 




V /dw 






(17) 

which plays the same role as Gc{K,lo) in the DCA loop. 
Once Gtyp is calculated from Eg. 1171 the self-consistency 
steps are the same as those in the multiband DCA de¬ 
scribed in the previous section: we calculate the coarse 
grained lattice Green function using Eq. [TOl and use it to 
update the hybridization function matrix of the effective 
medium via Eq. [TSl 

The proposed ansatz Eq. [11] has the following prop¬ 
erties. When the inter-band hopping and disorder 
potential 17“^ vanish, it reduces to single-band TMDCA 
for Ih independent bands, since all the off-diagonal el¬ 
ements of the Green functions vanish. When disorder 
is weak, all the 17““ are small so the distribution of the 
LDOS becomes Gaussian with equal linear and geometric 
average so it reduces to DCA for a multiband disordered 
system. 

When convergence is achieved, we use the total TDOS 
plyp{uj) to determine the mobility edge which is calcu¬ 
lated as the trace of the local TDOS matrix 


plllico) = Tr 


IT PtypjK,^) 


N, 


K 


Ptypi^)- ( 18 ) 

VCK —^ 


This construction of the order parameter may not seem 
very physical as the typical value of the LDOS should 
serve as the order parameteri^J^^, and the LDOS for the 
multiband system is the sum of the 4 bands in the 
local site basis P^^ Therefore, the 

real order parameter should be the typical value of 
defined as the geometric average of the total LDOS, 
exp(-^ log/5*°*) which is invariant under local uni¬ 
tary transformations and is not equal to the defined 
in Eq. [THl 

However, Eq. [TH] should also be a correct order param¬ 
eter as long as it vanishes simultaneously with the typ¬ 
ical value of and we show this in Appendix]^ By 


considering the distribution of the LDOS in each band. 
Appendix shows that when localized states mix with 
extended states the system is still extended, which is con¬ 
sistent with Mott’s insight about the mobility edge^. In¬ 
tuitively, this makes sense as when all the distributions 
of /9““ are critical then the typical values must behave 
as \V-Vc\^^ near the transition, and so their sum must 
as well. If one is not critical (on the metallic side) then 
Eq. [18] will not vanish as \V — Vc\^'', a.s expected. 

To test our multiband typical medium dynamical clus¬ 
ter approximation formulation, we apply it to the spe¬ 
cific case of a two band model, unless otherwise stated in 
Sec. mil Throughout the discussion of our results below, 
we denote a as a and /3 as b. 


III. RESULTS 


A. Two band model 


As a specific example, we test the generalized DCA 
and TMDCA algorithms for a three-dimensional system 
with two degenerate bands (ab) described by Eq. (Tj In 
this case, both the hopping and disorder potential are 2 
X 2 matrices in the band basis given by 

( ioa j.ab \ 

^ba ^bb j , ( 19 ) 


and 


fo* = 


^aa 
^ i 

yba 

I 



( 20 ) 


respectively. The intra-band hopping is set as = 
= 1, with finite inter-band hopping Here, the 
hopping matrix is defined as dimensionless so that the 
bare dispersion can be written as = tsk with Sk = 
—2t[cos{kx) + cos{ky) + cos{kz)] in three dimensions. We 
choose 4t = 1 to set the units of energy. We consider 
the two bands orthogonal to each other, where the local 
inter-band disorder 17“^ vanishes and the randomness 
comes from the local intra-band disorder potential I7““ 
that follow independent binary probability distribution 
functions with equal strength, I7““ = 17^^. Since the two 
bands are degenerate and the disorder strength for each 
band is also identical, the calculated average DOS will be 
the same for each band, so we only plot the quantities for 
one of the bands in the following results, as it is enough 
to characterize the properties of the system. 
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In our formalism, in order to disorder average instead 
of performing the very expensive enumeration of all dis¬ 
order configurations, which scales as 2^'^=, we perform 
a stochastic sampling of configurations which greatly re¬ 
duces the computational cost^. This is so we can study 
larger systems. For a typical Nc = 64 calculation, 500 
disorder configurations are enough to produce reliable re¬ 
sults and this number decreases with increasing cluster 
size. 

We first compare the ADOS and TDOS at various dis¬ 
order strengths with a hxed inter-band hop¬ 

ping = 0.3, for different cluster sizes Nc in Fig. [T] 
Our TMDCA scheme for Nc = 1 corresponds to the ana¬ 
log of the TMT for two-band systems, and the ADOS is 
calculated with the two-band DC A. To show the effects 
of non-local correlations introduced by finite clusters, we 
present data for both Nc = 1 and A"c > 1. We can clearly 
see that the TDOS, which can be viewed as the the or¬ 
der parameter of the Anderson localization transition, 
gets suppressed as the disorder increases . By comparing 
the width of the extended state region, where the TDOS 
is finite, we can see that single site TMT overestimates 
localization. 


aa bb , ab „ _ _ .ab „ 

Nc=l t =t =1, t =0.3, V =0 Nol 



FIG. 1. Evolution of the ADOS and TDOS at different dis¬ 
order strengths for Nc = 1 (left panel) and Nc > 1 

(right panel) for fixed t®*” = 0.3. For small disorder, the ADOS 
and TDOS are almost identical. The TDOS is suppressed as 
the disorder increases. The extended states region with finite 
TDOS for Nc = 1 is narrower than the results of Ac > 1 
which indicates that the single-site TMT overemphasizes lo¬ 
calization. 

From Fig. [U we see that the results of TMDCA for 
Nc = 64 and Nc = 216 are almost on top of each other, 
showing a quick convergence with the increase of cluster 
size. To see this more clearly, we plot in Fig. [2] the TDOS 
at the band center for two different disorder strengths 
and various cluster sizes. We see that the results for 
both cases converge quickly with cluster size. Faster 
convergence (around Nc = 38) is reached for the case 
further away from the critical region (t/““ = = 0.6) 


than for the one closer (I/““ = = 0.7) where con¬ 

vergence is reached around Nc = 98. This is expected 
due to the critical slowing down close to the transition. 
To further study the convergence, we also plot in Fig. [3] 
the TDOS at the band center as a function of disor¬ 
der strength (I/““ = for several Nc- The criti¬ 

cal disorder strength is dehned by the vanishing of the 
TDOS(a; = 0). The results show a systematic increase 
of the critical disorder strength as Nc increases, and the 
convergence is reached at Nc = 98 with the critical value 
of 0.74. 



FIG. 2. Evolution of the TDOS at the band center {u> = 0) 
with increasing cluster size for two different sets of parameters 
with r“ = = 1.0, = 0.3, F“'’ = 0.0, = 

0.6, 0.7. The former has faster convergence (around Nc = 38) 
than the latter (around Nc = 98), due to the critical slowing 
down closer to the transition region. 



with increasing cluster size, for = t*'*’ = 1.0, 1°"^ = 0.3, 

= 0.0. For Nc = 1, the critical disorder strength is 0.65 
and as Nc increases, it increases and converges to 0.74 for 
Nc = 98. 

To study the effect of inter-band hopping we calcu- 
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late the disorder-energy phase diagram for the case with 
vanishing and finite = 0.3 in Fig. S] The mobility 
edge is determined by the energy where the TDOS van¬ 
ishes. By comparing the left and right panels, we can see 
that introducing a finite makes the system more dif¬ 
ficult to localize, causing an upward shift of the mobility 
edge. The single site TMT overestimates the localized re¬ 
gion compared to finite cluster results. We also compare 
our results with those from the TMM to check the 
accuracy of the mobility edge calculated from TMDCA. 
For the TMM, the Schrddinger equation is written in 
terms of wavefunction amplitudes for adjacent layers in 
a quasi-one dimensional system, and the correlation (lo¬ 
calization) length is computed by accumulating the Lya¬ 
punov exponents of successive transfer matrix multiplica¬ 
tions that describe the propagation through the system. 
All TMM data is for a 3d system of length L = 10® 
and the Kramer-MacKinnon scaling parameter A{V,M) 
is computed for a given disorder strength V and “bar” 
width M. The transfer matrix is a 2Mlb x 2M/f, ma¬ 
trix. The system widths used were M = [4—12]. The 
critical point is found by identifying the crossing of the 
A{M)vs.V curves for different system sizes. The transfer 
matrix product is reorthogonalized after every five mul¬ 
tiplications. 



FIG. 4. Disorder-energy phase diagram for vanishing (left 
panel) and finite = 0.3 (right panel). We compare the 
mobility edge obtained from the TMT {Nc = 1), TMDCA 
(Ac = 64 and 216) and TMM. Parameters for the TMM data 
are given in the text (the TMM data for tab = 0.0 is repro¬ 
duced from -®). A finite increases the critical disorder 
strength, indicating that 1“*” results in a delocalizing effect. 
The single site TMT overestimates the localized region. 

To see the effect of inter-band hopping more directly, 
we now consider increasing while keeping the disor¬ 
der strength fixed = F®® = 0.71), and study the 
evolution of the mobility edge (Fig. [5]). The localized 
region around the band center starts to shrink as is 
increased, leading to a small dome-like shape with the 
top located at = 0.2. This shows that increasing 


delocalizes the system which is reasonable since increas¬ 
ing effectively increases the bare bandwidth. 



FIG. 5. Evolution of the mobility edge as increases, while 
F““ and F*'*’ are fixed. The results are calculated for Ac = 64. 
A dome-like shape shows up around the band center, signaling 
the closing of the TDOS gap. Parameters for the TMM data 
are given in the text. 

To further benchmark our algorithms, we plot the 
ADOS and TDOS calculated with two-band DCA and 
TMDCA together with those calculated by the KPM22ri22 
(Fig. [6]). In the KPM analysis, the LDOS is expanded 
by a series of Chebyshev polynomials, so that the ADOS 
and TDOS can be evaluated. The details for the imple¬ 
mentation of KPM are well discussed in Ref. HO and the 
parameters used in the KPM calculations are listed in 
the caption of Fig. [6l The Jackson kernel is used in the 
calculations^. As shown in the plots, the results from the 
generalized DCA and TMDCA match nicely with those 
calculated from the KPM. 

The excellent agreement of the TMDCA results with 
those from more conventional numerical methods, like 
KPM and TMM, suggest that the method may be used 
for the accurate study of real materials. 


B. Application to KyFe2-a:Se2 

Next, we demonstrate the method with a case 
study of Fe vacancies in the Fe-based superconductor 
K 3 ,Fe 2 -ySe 2 , which has been studied intensely because 
of its peculiar electronic and structural properties. Early 
on it was found that there is a strong y/b x Vb or¬ 
dering of Fe vacancies^. Later it was discovered that 
this material also contains a second phas o®^’®® . It is 
commonly speculated that the second phase is the one 
that hosts the superconducting state and the phase with 
the X vacancy ordering is an antiferromagnetic 
(AFM) insulator. Recent measurements of the local 
chemical compositio n®®'®^ have determined that the sec¬ 
ond phase also contains a large concentration of Fe va¬ 
cancies (up to 12.5%). However, these Fe vacancies are 
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FIG. 6. Comparison of ADOS and TDOS calculated with 
DCA, TMDCA and KPM with fixed disorder strength 1/““ = 
V'’*’ = 0.8 and various values of inter-band hopping 1“*". The 
KPM uses 2048 moments on a cubic lattice of size 48® and 
200 independent realizations generated with 32 sites randomly 
sampled from each realization. 

not well ordered since no strong reconstruction of the 
Fermi surface^^— was observed by angle-resolved pho¬ 
toelectron spectroscopy (ARPES) experiments^!^. 

Interestingly, with such a disordered structure, this 
material hosts a relatively high superconducting transi¬ 
tion temperature of 31 K at ambient pressure^. It was 
the first Fe-based superconductor that was shown from 
ARPES^i^ to have a Fermi surface with electron pock¬ 
ets only and no hole pockets, apparently disfavoring the 
widely discussed pairing symmetry^ in the Fe-based 
superconductors. K 2 ,Fe 2 -ySe 2 is also the only Fe-based 
superconductor whose parent compound (with perfectly 
ordered Fe vacancy) is an AFM insulator— rather than 
a AFM bad metal. Furthermore from neutron scatter¬ 
ing^, it has been observed that the anti-ferromagnetism 
has a novel block type structure with a record high Neel 
temperature of T^v = 559K and magnetic moment of 
S.Sl/rs/Fe. Such a special magnetic structure is obvi¬ 
ously not driven from the nesting of the simple Fermi 
surface, but requires the interplay between local moments 
and itinerant carriers present in the normal state^i^. 

Given that Fe vacancies are about the strongest possi¬ 
ble type of disorder that can exist in Fe-based supercon¬ 
ductors and given that the Fe-based superconductors are 
quasi two-dimensional materials, it is natural to specu¬ 
late how close the second phase is to an Anderson insu¬ 
lator. If it is indeed close, this would have interesting 
implications for the strong correlation physics and the 
non-conventional superconductivity in these compounds. 

To investigate the possibility of Anderson localiza¬ 
tion in the second phase of Ka;Fe 2 -ySe 2 we will employ 
TMDCA on a realistic first principles model. To this end 
we use Density Functional Theory (DFT) in combination 
with the projected Wannier function technique^ to ex¬ 
tract the low energy effective Hamiltonian of the Fe-d de¬ 


grees of freedom. Specifically we applied the WIEN2K^ 
implementation of the full potential linearized augmented 
plane wave method in the local density approximation. 
The k-point mesh was taken to be 10 x 10 x 10 and the 
basis set size was determined by RKmax=7. The lattice 
parameters of the primitive unit cell (c.f. Fig. [3(b)) are 
taken from Ref. l33l The subsequent Wannier transfor¬ 
mation was defined by projecting the Fe-d characters on 
the low energy bands within the interval [-3,2] eV. For 
numerical convenience, we use the conventional unit cell 
shown in Fig. IZKa) which contains 4 Fe atoms. Since 
there are 5 d orbitals per Fe atom, we are dealing with a 
20-band problem. To simulate the effect of Fe vacancies 
we add a local binary disorder with strength V and Fe 
vacancy concentration Ca- 

pm = Ca6{V^ - E) + (1 - Ca)<5(E)- (21) 

We set the disorder strength to be E = 20eE, much 
larger than the Fe-d bandwidth, such that it effectively 
removes the corresponding Fe-d orbitals from the low en¬ 
ergy Hilbert space. This will capture the most dominant 
effect of the Fe vacancies. The Fe concentration is taken 
to be Cq = 12.5%, which is the maximum value found in 
the experiments. 


(a) conventional unit cell (b) primitive unit cell 


FIG. 7. Grystal structure of KFe 2 Se 2 . 

Fig. m presents the ADOS and TDOS, obtained from 
our multiband TMDCA for which we considered two 
cluster sizes Nc = 1 and Nc = 2y/2 x 2-\/2 x 2 = 16. 
Consistent with the model calculations presented in the 
previous sections, we find that the TMT (W = 1) 
tends to overestimate the localization effects compared to 
TMDCA results (W = 16). While the TMT shows local¬ 
ized states within [0.6,1.1] eV, the TMDCA for W = 16 
finds localized states in the much smaller energy re¬ 
gion [1.0,1.1] eV instead. Apparently a concentration of 
Ca = 12.5% is still too small to cause any significant lo¬ 
calization effects despite the strong impurity potentials 
of the Fe vacancies and the material being quasi-two di¬ 
mensional. To determine the chemical potential we con- 


























sider two fillings. The first filling of 6.0 electrons per 
Fe corresponds to the reported K 2 Fe 7 Se 8 phased. Since 
strong electron doping has been found in ARPES exper- 
iments^i^, we also consider a filling of 6.5 electrons per 
Fe. The latter would correspond to the extreme case of 
no vacancies. Clearly for both fillings the chemical po¬ 
tential remains energetically very far from the mobility 
edge, and thus far from Anderson insulating. 



FIG. 8. The average and typical density of states of KFe 2 Se 2 
with 12.5% Fe vacancy concentration calculated by multiband 
DCA and TMDCA with cluster size Ac = 1 and Ac = 16, 
compared with the average density of states of the clean (no 
vacancy) KFe 2 Se 2 . 


IV. CONCLUSION 

We extend the single-band TMDCA to multiband sys¬ 
tems and study electron localization for a two-band 
model with various hopping and disorder parameters. We 
benchmark our method by comparing our results with 
those from other numerical methods (TMM and KPM) 
and find good agreement. We find that the inter-band 
hopping leads to a delocalization effect, since it gradu¬ 
ally closes the a; = 0 disorder induced gap on the TDOS. 
A direct application of our extended TMDCA could be 
done for disordered systems with strong spin-orbital cou¬ 
pling. Combined with electronic structure calculations, 
our method can be used to study the electron localiza¬ 
tion phenomenon in real materials. To show this, we 
apply this approach to the iron selenide superconductors 
K 3 ,Fe 2 _ySe 2 with Fe vacancies. By calculating the TDOS 
around the chemical potential, we conclude that the in¬ 
sulating behavior of its normal state is unlikely due to 
Anderson localization. This method also has the abil¬ 
ity to include interactions^^, and future work will involve 
real material calculations that fully treat both disorder 
and interactions. 

Acknowledgments- This work is supported in part by 
the National Science Foundation under the NSF EP- 


SCoR Cooperative Agreement No. EPS-1003897 with 
additional support from the Louisiana Board of Re¬ 
gents (YZ, HT, CM, CE, KT, JM, and MJ). Work 
by TB was performed at the Center for Nanophase 
Materials Sciences, a DOE Office of Science user fa¬ 
cility. This manuscript has been authored by UT- 
Battelle, LLC under Contract No. DE-AC05-00OR22725 
with the U.S. Department of Energy. WK acknowl¬ 
edges support from U.S. Department of Energy, Of¬ 
fice of Basic Energy Science, Contract No. DEAC02- 
98CH10886. The United States Government retains and 
the publisher, by accepting the article for publication, 
acknowledges that the United States Government re¬ 
tains a non-exclusive, paid-up, irrevocable, world-wide 
license to publish or reproduce the published form of this 
manuscript, or allow others to do so, for United States 
Government purposes. The Department of Energy will 
provide public access to these results of federally spon¬ 
sored research in accordance with the DOE Public Access 
Plan (http://energy.gov/downloads/doepublic- access- 
plan). This work used the high performance computa¬ 
tional resources provided by the Louisiana Optical Net¬ 
work Initiative (http://www.loni.org), and HPG@LSU 
computing. 

Appendix A: The order parameter defined in Eq. 1181 

We know the system is localized if the distribution of 
the total LDOS is critical, having a probability distribu¬ 
tion p(/o““ -I- which is highly skewed with a typical 
value close to zero. So if we can show that this is true if 
and only if both and pf’ are critical, then the criti¬ 
cal behavior is basis independent and we can choose any 
particular basis and use the order parameter defined by 
Eq. [18] to study the localization transition. 

To show this is true, we consider two probability dis¬ 
tribution functions pi{xi) and ^ 2 ( 2 ^ 2 ). The probability 
distribution function for A = a;i -I- 0:2 is 

P{X) = f pi{x)p 2 {X - x)dx, (Al) 

Jo 

and we want to show P{X) is critical if and only if both 
Pi(a;i) andp2(a:2) are critical. 

1. Sufficiency 

If both pi{x) and ^ 2 ( 2 ;) are critical, then both pi{x) 
and P 2 {x) are dominated by the region 0 < x < S where 
(5 —>■ 0+. The contribution to the integral in P(A) mainly 
comes from the region 0 < a: < d and 0 < A — a; < <5 
which is niax{X — 6,0) < x < min{5,X). Since 5 is 
inhnitesimal, we can assume X > 5, and then we have 
X — 6 < X < 6. To maximize P(X), we want this region 
to be as big as possible, so we want <5 — (A — 6) = 26 — X 
to be as big as possible which means A must be smaller 











2. Necessity 
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than 2S —>■ O'*". Thus, P{X) is also critical with the 
typical value around 25 which is infinitesimal. 


We now consider the case where one of the distribu¬ 
tions is not critical. Without loss of generality, we assume 
P 2 {x) is not critical and is peaked at some finite value xq- 
We calculate 


P{xo)-P{5) 


nXQ nS 

/ pi{x)p 2 {xq — x)dx — / pi{x)p 2{5 — x)dx 
Jo Jo 

rS /-xo 

/ pi{x)[p 2 {xq - x) - P 2 {S - x)]dx + pi{x)p 2 {xo - x)dj 
Jo Js 


(A2) 


The first term is positive since p 2 {x) is peaked around xq 
and 5 Xo- The second term is positive obviously, so 
P(xo) > P(5). Therefore, P(X) is not critical. 

In this way we argue that P(X) is critical if and only if 
both pi{xi) and P 2 {x 2 ) are critical. In other words, when 


the localized states hybridize with extended states, only 
extended states remain which is exactly Mott’s insight 
about the mobility edge^i. The generalization to the 
multiple band case is trivial. 
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